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The low energy n — 3 H elastic cross sections near the resonance peak are calculated by solving the 
4-nucleon problem with realistic NN interactions. Three different methods - Alt, Grassberger and 
Shandas (AGS), Hyperspherical Harmonics and Faddeev-Yakubovsky - have been used and their 
respective results are compared. We conclude on a failure of the existing NN forces to reproduce 
f"**^ ' the n — 3 H total cross section. 

CN ' PACS numbers: 21.45.+V, 11.80.Jy, 25.40.-h, 25.10.+S 
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I. INTRODUCTION 

The four nucleon (4N) system represents a qualitative jump in complexity relative to the A — 3 case, as has been 
put forward in several papers 0,0 11 011110 HOI- It becomes already obvious when comparing the experimental 
n — p, n — d and n — 3 H cross sections displayed in Fig. Q ; the smooth behaviour of A = 2 and A = 3 curves contrasts 
with the non trivial structure manifested in A — 4. This structure, in the energy region around E cm «3 MeV, is 
commonly associated with negative parity, isospin T = 1 resonance states |l0| . Actually, A — 4 is the smallest nuclear 
— ■«. ■ system which exhibit characteristic nuclear properties - such as saturation in the binding energy - and in which the 
\ simplest nuclear (n + 3 He — * p + 3 H) and fusion (d + d — ► n + 3 He) reactions take place. Moreover, many reactions 
involving four nucleons, like p + 3 He — > 4 Hc + v c + c + (the hep process), are of extreme astrophysical interest, as they 
play important roles in solar models and big-bang nucleosynthesis; the hep process, for instance, is the source of the 
^ . highest energy neutrinos from the Sun. 

The study of the 4N system is particularly interesting as a "theoretical laboratory" to test new models of the nuclear 
J> , force. Unlike the A = 3 systems, A = 4 shows a delicate and rich structure of excited states in the continuum E3 
whose position and width depends critically on the underlying nucleon- nucleon (NN) interaction In fact, the 
effect of (i) the NN P-wave and of (ii) the three-nucleon (3N) force are believed to be larger than in the A = 2 or 3 
systems. Moreover, it is the simplest system where the 3N interaction in the channels of total isospin T = 3/2 can be 
studied. Therefore it is of the utmost importance to have reliable few-body techniques powerful enough to deal with 
this system. Ever increasing computer power, development of novel numerical methods, and significant refinements 
of well-established techniques have come together to show that the solution of the four-nucleon bound state problem 
with realistic Hamiltonians is reliable. 

In Ref. [HI the binding energies and other properties of the a-particle were studied using the AV8' E2 NN 
interaction; several different techniques El Q [ll [ll El El El HI IU IH |H HHI produce results in very close 
agr eement with each other (at the level of less than 1%). More recently the Faddeev-Yakubovsky (FY) equations [51 
1261 ] were solved for different realistic NN + 3N force models. These calculations are fully converged in terms of NN 
partial waves and the results agree with identical works using the Green's function Monte Carlo (GFMC) |2l| or 
Hypherspherical Harmonics (HH) methods |Tsj |. 

In this paper, we consider the problem of solving the 4./V scattering problem with realistic Hamiltonians. As 
discussed before, the A = 4 system is very rich and several elastic and rearrangement channels are indeed possible 
already at low energy, depending on the total charge (n — 3 H, p — 3 H <-> n — 3 He ^> d — d, p — 3 He). Among them, the 
n — 3 H is the simplest one: free from Coulomb interaction and, to a very good approximation, a pure isospin T = 1 
state. However, it presents a rich dynamics since a resonance structure is visible at center of mass energy E « 3 MeV, 
as already shown in Fig. QJ. 

Up to now only two of the techniques used in bound state calculations (FY and HH) as well as methods based on 
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the resonating group model (RGM) [27| and on the solution of Alt, Grassberger and Sandhas (AGS) equations l^.[23l 
have been employed to attack this problem with realistic Hamiltonians (for a review of earlier results, see Ref. jlCj')- 
Some previous FY calculations in configuration space suggested that realistic NN interactions could fail in reproducing 
the experimental total and differential cross sections 0,11, E| . These calculations were however limited to a relatively 
small number of partial waves in the expansion of FY amplitudes, thereby such lack of convergence could be advocated 
to explain the failure. Using AGS equations, the author of Ref. was able to substantially increase the number 
of partial wave amplitudes and obtain a fairly good agreement for the total cross section at the resonance peak. 
However, these calculations were based on a rank-one separable expansion of the 2-body t-matrix, which also raises 
some doubts on the reliability of the results. The calculations performed by expanding the scattering wave functions on 
the correlated Hyperspherical Harmonic (CHH) basis were able to obtain a reasonable estimate of the various phase- 
shift and mixing-angle parameters However, the problem of convergence could not be solved completely due to 
numerical difficulties when adding more CHH components. For this reason, a new expansion on the (uncorrelated) 
HH basis is considered here. 

The results obtained using these techniques were so far at rather large variance between each other. Clearly, 
this situation should be clarified before questioning the ability of present NN + 3N force models to describe the 
experimental data beyond the binding energy of 4 He, for which they have proven to be rather successful plj . This is 
the purpose of the present paper in which we compare low energy n — 3 H scattering results obtained by three different 
groups, using independent methods to solve the four-body equations. They concern AGS calculations Q based on 
rank-one expansion of the NN t-matrix, variational solutions of the Schrodinger equation using HH method 0] and 
the solutions of FY equations [2j| [3(| with a substantial improvement in the number of partial waves j^l); AV18 [12 
and NIJM-II 33] NN interaction models are used in present calculations. 

We have considered the center of mass energies E = 0.40, 0.75, 1.50, 2.625, 3.0 MeV. They constitute a good set of 
data for interpolating other values of the phase shifts and allow to determine the low energy parameters (scattering 
length and effective range). In addition, values E = 0.75,1.50,2.625 MeV correspond to the measured differential 
cross sections. 

To describe scattering solutions, we will expand the different J 71 ' n— 3 H states in terms of the asymptotic hamiltonian 
channels \L, S; J v >, where L denotes the n — 3 H relative orbital angular momentum and S = s n + s t the total spin. 
The following J 17 states will be considered: 

|0,0;0+) 
c 1+ |0,l;l+)+rf 1+ |2,l;l+) 

|1,0;0-) 
Cl - |l,l;l-)+d 1 -|l,0;l-) 



|0+) = 
|l + > = 

|o-) = 

m = 
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|2") = c 2 - |l,l;2-)+d2-|3,l;2-)| (1) 

Coefficients c and d are fixed by the dynamics and make the solutions eigenstates of the S- matrix. 

As it will be demonstrated later, these states are the only relevant ones for the low energy scattering considered 
here. Moreover, the coupling L L + 2 is very weak and is neglected for J n = 2~ state. 

II. METHODS 

We briefly describe the methods used for solving the 4N problem. Two of them (FY and HH) work in configuration 
space while AGS in momentum space. 

A. AGS equations 

The starting point involves the AGS equations [2^] for the transition operators involving all (2N)+(2N) and N+(3N) 
channels. For local NN potentials such equations are three-vector variable integral equations which after partial wave 
decomposition reduce to a set of coupled equations in three continuous scalar variables. Since scattering calculations 
require a great number of channels for convergence, we follow an approach based on the separable representation of 
subsystem amplitudes in order to reduce the equations to two or one continuous variable. The integral equations we 
use are the same as in Ref. [34| and result from the modified AGS equations j3f| after one has: (a) represented the 
original NN t-matrix by an operator of rank one; (b) represented the resulting 3N t-matrix by a finite rank operator 
and taken as many terms as needed for convergence. Since in the modified AGS equations the 2N+2N subamplitudes 
are expressed in terms of a convolution integral involving two non-interacting pair-propagators, as first proposed in 
Ref. |36( , the sole approximation in this approach involves a rank one representation of the 2N t-matrix which may be 
obtained from the well-known method of Ernst, Shakin, and Thaler (EST) [33. The multi-term representation of the 
3N t-matrix is done using the Energy Dependent Pole Expansion (EDPE) method developed in Ref. (3^ . This latter 
representation for the 3N t-matrix is well under control since one may check the convergence rate of 4N observables 
for increasing rank in the 3N t-matrix. 

This method was first used in Ref. [23| to calculate the binding energy of 4 He and later confirmed to be accurate 
by the exact work of Kamada and Glockle 0|. More recently the results of our calculations for n — 3 H elastic 
scattering were shown to agree with the exact results of the Grenoble group Q, for both Malfliet-Tjon (MT) and 
AV14 potentials taken in 2N partial waves with total angular momentum j < 1 + ( 1 So, 3 Si — 3 Di). In the present 
calculations we extend the work in Ref. [|| by increasing the number of 2N partial waves to j < 2. Therefore we add 
3 F 2 , l D 2 and 3 D 2 to the partial waves 1 S , 3 Si - 3 Di, 1 P 1 , 3 P , 3 Pi and 3 P 2 already included in Ref. @. Moreover 
we limit the particle-pair orbital angular momentum £ y < 2 and include all J3 and £ z < 2 for a given total four-body 
angular momentum J; J 3 is total angular momentum of any given three-body subsystem and £ z is the orbital angular 
momentum of that subsystem relative to the fourth particle. The results we find here are consistent with those 
obtained in Ref. and converged relative to the number of terms in the finite rank expansion of the underlying 3N 
t-matrix. The number of terms in the expansion range from four to six depending on J3. 

B. Faddeev-Yakubovsky method 

In the case of four identical fermions, interacting via a pair- wise potential V, the FY equations result into a set of 
two integrodifferential equations, coupling two FY components, namely K and H: 

{E-H Q -V)K = V{P+ +P-)[{1 + Q)K + H] (2) 
(E-H -V)H = VP[{l + Q)K + H] (3) 

with P + , P~ , P and Q being particle permutation operators: 

P + = (P-)- = P23-P12; Q = P = PiiP2i = P24P13, (4) 

and e is a Pauli factor for exchange of two identical particles, which in case of fermions is -1. The wavefunction is 
given by: 

* = [1 + (1 + P + + P-)Q] (1 + P + + P-)K + (1 + P + + P~)(l + P)H (5) 
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Each FY component F~(K,H) is considered as a function, described in its proper set of Jacobi coordinates x, y, z, 
denned respectively by 



xk = r 2 -n _ _ 

Vk = V 3 ^ 2^) VH = _ n-r 3 (6) 



2 V' 4 3 

and expanded in angular variables for each coordinate according to 



(xyz\F) = Y a (x,y,z). (7) 

* — ' xyz 



The quantities F a are called regularized FY amplitudes and Y a are tripolar harmonics, containing spin, isospin and 
angular momentum variables. The label a holds for the set of 10 intermediate quantum numbers describing a J 71 ", T = 1 
state . They are defined in a j-j coupling scheme as 



K = {[{t^rMrJ^T^ ® \ ( l x(siS 2 )aJ jx (lyS 3 ) jy ^^(l z S4,) jz 



' 3y 



where and tj are the spin and the isospin of the individual particles and T, J the isospin and total angular 
momentum of the four-body system. Each of the N c = Nk + Nh amplitudes in the expansion J7J) is labelled by 
12 quantum numbers, which are further conditioned by the antisymmetry properties (— ) <T =+ r =+^ — £ for K and 

The boundary conditions for the 1+3 scattering problem are implemented by imposing at large enough value of z 
the Dirichlet-type condition 

K(x,y,z) = t(x,y) 
H(x,y,z) = 

t(x, y) being the triton Faddeev component with quantum numbers (^(siS2) CTx ),■ (lySa)j y ■ They ensure a solution 
which, e.g. for a relative n — 3 H S-wave, behaves asymptotically like 



where S is the 



K(x,y,z) ~ t(x, y) sin (gz + 5) 
3 H phase shift and q, the conjugate momentum of the z- Jacobi coordinate in K-amplitudes, is 



related to the center of mass n — H kinetic energy E cm and the physical momentum k by 



7-Elab — — q 

4 771 



ft? n 2ft 2 , 



3 777 



(8) 



The convergence is reached in respect to expansion J7J). 



FY calculations have been performed in the j — j coupling scheme. The following truncations in the partial wave 
expansion of amplitudes were used: (i) Vnn waves limited to l x < 3, always including tensor-coupled partners, i.e. 
involving the set ^o, 3 SD X , 1 P 1 ? P , 3 PF 2 , 3 i\, l D 2 , 3 DG 3 , 3 D 2 , X F 3 , 3 FH 4 3 F 3 and (ii) l x + l y + l z < 10. 

The convergence was studied as a function of j yz =max.{j y ,j z ) for K-like components and j yz =max(jy,l z ) for H-like, 
starting with j yz = 1. In the second column of the Table[I]we present the values of the relevant phase shifts for AV18 
model at the peak energy E cm =3.0 MeV. The couplings L «-> L + 2 in the asymptote of the 1 + (L = <-> L = 2) 
and 2 _ (L = 1 <-> L = 3) states are found to be very small (see next section) and have been omitted here. For the 
1~ state, we give the eigenphaseshifts - Si(l~), 5 2 (1~) - and the mixing parameter e. One can remark that results 
displayed on Table [I] converge pretty well as a function of j yz . For the most of the phaseshifts three-digit accuracy 
is reached already with j yz = 2. The most unstable were found J T = 2~ phascshift and J w = 1~ mixing parameter, 
which were requiring inclusion of j yz — 3 amplitudes to converge at a one percent level. 
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C. Hyperspherical harmonic method 



In the HH method, the wave function is written as the sum of an asymptotic part, which form is known (except 
for the parameter of the S- matrix), and an internal part, which is expanded in HH functions. Explicitly, the wave 
function ^x +3 J " describing an — 3 H scattering state with incoming orbital angular momentum L and channel spin S 
(S = 0, 1) coupled to total angular JJ Z , is expressed as 

where vanishes in the limit of large intercluster separations, and hence describes the system in the region where 

the particles are close to each other and their mutual interactions are strong. On the other hand, ^! 1 ^ SJJz describes 
the relative motion of the two clusters in the asymptotic regions, where the n — 3 H interaction is negligible and can 
be decomposed as a linear combination of the following functions 

^ 4 

nisjj. = "/J E MjIp)]s ® Y L (*ij\ Jh (-fRirMkn) ± ij L (kn)) , (10) 

where r$ is the distance between the neutron (particle i) and 3 H (particles jlp), k is the magnitude of the relative 
momentum between the two clusters defined in Eq. (HJ , 03 is the 3 H wave function, which has been obtained by solving 
the corresponding three-body problem using the pair-correlated HH expansion |39| and jh an d TJl are the spherical 
Bessel functions of the first and second kind, respectively. The function fn(ri) has been introduced to regularize 
VL{kri) at small r,, and /ii(ri) — > 1 as r% is large, thus not affecting the asymptotic behavior of ^if^ Jz . Note that 
for large values of fcfj, 



cos(fcr, - Ltt/2) ± i sin(fcri - Lir/2) ex P 



±i(kn - Ltt/2) 



-fn(ri)yL(kri) ± y'i(feri) -> ; = - ; - . (11) 



kri kri 

L LSJJz y^LSJJz 



Therefore, ^Js/j z {^lsjj^ describes in the asymptotic regions an outgoing (ingoing) n — 3 H relative motion. Finally, 



(12) 



SlL>SsS'^LSJJ„ - S LS,L'S'( k )^LSJJ, 

US' L 

where the parameters S"[ s L , s , (k) are the 5-matrix elements which determine phase-shifts and (for coupled channels) 
mixing angles at the energy (2/3)fc 2 /m . Of course, the sum over L' and 5" is over all values compatible with a given 
J and parity. 

The "core" wave function ^ J c Jz ^ has been here expanded in the HH basis as for the bound-state wave function [is| . 

N c n 2 +n 3 <N(a) 

^ J C' 7 ' = X! ^2 u an 2 n 3 {p)y a n 2 n 3 (xK,yK,Z K ) , (13) 

a— 1 n2.n 3 — 



where p is the hyperradius, p = y/ x\ + y\ + z 2 K . The known functions y a n 2 n 3 {xK,yK,ZK) are given as the anti- 
symmetrized product of spin-isospin states and HH functions |40|, the latter being given by the product of tripolar 
spherical harmonics and two Jacobi polynomials of indexes ri2 and given, respectively, in terms of "hyperangles" 
defined by 

COS <p 2 = cos ip 3 = . (14) 

V z k + Vk P 

The channels index a denotes collectively the orbital angular momenta, spin and isospins quantum numbers. For a 
given channel a, the HH states with indexes n-i and 77,3 chosen in the range < n% + TI3 < N(a) are included in 
the expansion. The values of the number of channels iV c and the various N(a) are increased until the desired degree 
of convergence in the quantity of interest is obtained (for the choice of which channels to include see the discussion 
below) . 

The functions u arL2n3 (p) occurring in the expansion of tyg Jzn and the matrix elements S£ s L , s , (k) are determined 
by making the functional 

&s,L>s>(k)] = SisMS'ik) - ^<<f JJ *|ff - E, - T cm \*^i J *) (15) 
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stationary with respect to variations in the S^ s L , s , and u ari2n3 (Kohn variational principle). Here £3 is the 3 H ground- 
state energy. By applying this principle, a set of second order differential equations for the functions u an2n3 (p) are 
obtained. By replacing the derivatives with finite differences, a linear system is obtained, which have been solved using 
the Lanczos algorithm. The procedure is very similar to that outlined in the appendix of Ref. |4l| and it will not be 
repeated here. A large number of equation can be quite easily be solved. 

Let us now discuss briefly the choice of the HH states to be included in the expansion. As for the a-particle 
bound state [l8|, the "brute force" inclusion of HH states is not possible and one has to select the appropriate 
subset of functions. We have found convenient |l8| to separate the HH functions in classes depending on the value of 
C = £ x + £ y + £ z : channels with low values as possible of C have to be included in the expansion first. For each group 
of channels of a given £, the corresponding values of N(a) are increased until convergence is reached. Then, channels 
with a larger value of C are included in the expansion and so on. The advantage of this procedure is that the classes 
of channels with large values of C have small contributions. Consequently, the corresponding values of N(a) can be 
also taken smaller and smaller. 

An example of convergence is presented in Tabic [IJ in the last three columns. There, the phase shifts at E cm = 3 
MeV obtained by including in the HH expansion the channels with increasing values of C have been reported. For 
S-wave phase shifts, the inclusion of the C = channels is sufficient to have already rather good estimates. The 
contribution of the C > 4 channels is found to be negligible. For the P-wave phase shifts, the channels with C = 1 give 
the dominant contribution at low energies (i? c . m . < 1.5 MeV), but in the peak region they are clearly insufficient to 
obtain good estimates. Including the channels with C = 3, the HH predictions already differ less of 4% with respect to 
the final FY results. There exist a very large number of £ = 5 channels and we have limited ourselves to the inclusion 
of those (usually) most important, namely, those having t x — t y = 2 and l z = 1, plus permutations (this is sufficient 
to double the number of channels). However, the changes of the P-wave phase shifts are rather tiny. Therefore, we 
expect that the contribution of the disre gard ed channels C = 5 (and those with L > 5) to be practically negligible. 
This issue is currently being investigated |42l |. 

It must be noticed that in the present HH calculation only states constructed in terms of the Jacobi vectors of 
the set K have been considered (hereafter referred to set-K HH functions). Obviously, the HH functions of a given 
set of Jacobi vectors form a complete basis by themselves. However, the HH functions constructed with the set H of 
Jacobi vectors (set-H HH functions) should be more suitable for taking into account contributions from the {2 + 2} 
cluster structures. It is rather obvious that the inclusion of HH functions of both sets should speed up the convergence 
in constructing the state of the system. In the present calculation, however, we have included only the set-K HH 
functions since we were able to include a number of states sufficient to reach, in the studied cases, a satisfactory 
degree of convergence. 

TABLE I: Convergence of n — 3 H phase shifts and mixing parameter for AV18 potential at E cm =3 MeV. The coupling L — » L + 2 
in the asymptotics for the J n = 1 + and 2~ states have been neglected. See the text for explanations. 





FY 


HH 


r 


JVZ < 1 jyz < 2 jyz < 3 jyz < 4 


£=0£=2£=4 


0+ 


-70.11 -69.98 -70.00 -70.00 


-70.1 -69.8 -69.8 


it 


-63.61 -61.94 -61.97 


-62.6 -62.2 -62.2 


r 


jyz < 1 jyz < 2 jy Z < 3 jyz < 4 


C = l £ = 3 £ = 5 


(T 


26.28 23.86 23.60 23.58 


16.9 23.0 23.4 


ir 
ij 

e 


25.32 22.40 22.27 
18.00 40.70 40.88 
9.04 -43.52 -44.39 


20.7 22.0 22.2 
31.5 41.5 42.0 
-40.9 -44.9 -45.3 


2" 


43.82 45.41 45.44 


27.9 44.0 44.5 
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III. RESULTS 

The n — 3 H scattering lengths have already been presented in Ref. 0, |(| |4^. For completeness, in Table ITT1 we 
summarize the latest obtained values with the improved precision. 

In Table IIIII wc present phase shifts, mixing parameters and the corresponding total cross sections for the most 
relevant partial waves at some selected energies. The coupling L «-> L + 2 between the different asymptotics in the 
J n = 1 + , 2~ states turned out to be very small. In AGS and FY calculations they still have been taken into account 
for J 71 " = 1 + but neglected for J 71 — 2~ state. In fact, even at E cm = 3 MeV, the largest energy we have considered, 
L = 2 phaseshift in J 71 " = 1 + is only ~ 1° and contributes by only 0.01% in total cross section. Contribution of L = 3 
phaseshifts, which were neglected in J 77 = 2" calculations, is at least by one more order of magnitude smaller. 

TABLE II: Singlet ai(J 7r =0 + ) and triplet a 3 (J 7T = l + ) scattering lengths calculated with Nijm-II and AV18 NN interactions. 



Pot. 


ai (fm) 


a 3 (fm) 


a(0) (b) 


Method 




4.09 


3.70 


1.82 


AGS 


AV18 


4.25 


3.74 


1.88 


FY 




4.28 


3.73 


1.89 


HH 


Nijm-II 


4.09 


3.70 


1.82 


AGS 




4.25 


3.74 


1.88 


FY 



One can see that positive parity states are under full control by the different methods. Agreement between HH and 
FY results is perfect for positive parity phaseshifts (in the most of cases all three significant digits coincide, whereas 
the largest discrepancy does not exceed 0.3%). AGS results also agree by a few percent. 

The situation is more delicate for negative parity states, where resonances are present and various techniques suffer 
from slower convergence. Nevertheless FY and HH results stay in close agreement to better than few percent for 
phaseshifts, though it could be perfect if one or more convergence step is effected (see Table QJ. Still these small, 
convergence related, correlations do not shelter the excellent agreement in total cross sections. On the other hand 
disagreement with AGS results becomes important, reaching 30% for the phaseshifts in J" = 0~ state. All the 
negative parity phases are strongly overestimated by this method, thus resulting a larger total cross sections in the 
resonance region. 

We believe that the underlining reason for the disagreement with AGS results is due only to the rank-one separable 
expansion for 2-body t-matrix used in this calculations. This lowest order expansion may be enough to describe 
systems dominated by spherical symmetric wave functions (i.e. positive parity states) but fails for more complicated 
structures, not being able to account for strong compensations present in NN P-waves |43l |. which turn out to be 
sizeable in negative parity states CI- 
TABLE III: AV18 n - 3 H S- and P-wave phase-shifts (degrees) as a function of the cm. kinetic energy E cm (MeV). Numbers 
in brackets are the mixing parameters in degrees for corresponding scattering states. 





0+ 


It 


1 2 + 




0" 




2" 


<r(b) 




0.40 


-27.8 


-25.0 


-0.018 




2.24 


2.87 3.89 




3.56 


1.73 


AGS 




-28.8 


-25.1 


-0.0165 (0.199) 


1.81 


2.66 3.56 


(-39.6) 


3.38 


1.75 


FY 




-28.8 


-25.1 






1.78 


2.70 3.76 


(-44.2) 


3.32 


1.76 


HH 


0.75 


-37.3 


-33.6 


-0.072 




5.59 


6.69 9.95 




9.15 


1.79 


AGS 




-38.7 


-33.8 


-0.059 


(0.358) 


4.41 


6.11 8.85 


(-41.3) 


8.57 


1.78 


FY 




-38.7 


-33.8 






4.36 


6.08 9.32 


(-45.4) 


8.53 


1.79 


HH 


1.50 


-50.7 


-46.1 


-0.274 




14.5 


15.1 25.9 




24.3 


2.22 


AGS 




-52.8 


-46.2 


-0.330 


(0.670) 


11.5 


13.4 21.9 


(-42.9) 


22.5 


2.06 


FY 




-52.6 


-46.3 






10.8 


13.3 23.0 


(-45.6) 


22.1 


2.06 


HH 


2.625 


-63.7 


-58.3 


-0.649 




27.5 


24.7 44.8 




44.2 


2.51 


AGS 




-66.5 


-58.5 


-0.850 


(1.08) 


20.9 


20.7 37.3 


(-43.5) 


41.0 


2.24 


FY 




-66.3 


-58.7 






20.6 


20.5 38.6 


(-45.5) 


40.1 


2.24 


HH 


3.0 


-67.4 


-61.9 


-0.764 




30.9 


26.8 48.5 




48.5 


2.48 


AGS 




-70.0 


-62.0 


-1.180 


(1.21) 


23.9 


22.4 40.7 


(-43.5) 


45.4 


2.21 


FY 




-69.8 


-62.2 






23.4 


22.2 42.0 


(-45.3) 


44.5 


2.21 


HH 



The total cross sections are plotted in Fig. [21 and compared to the experimental values taken from Ref. [45j . The 
results of HH and FY (dashed line) are in agreement for three significant digits and thus are not distinguishable by eye. 



8 



The AGS cross sections (solid line) are slightly smaller at very low energies, whereas they provide considerably larger 
cross sections in the resonance region. The disagreement of theoretical results with experimental data at zero energy 
can be attributed to the 3 H underbinding by local NN interaction models. It can be improved by including a three 
nucleon force (3NF), i.e. the Urbana version IX (UIX), and thus reproducing the 3N binding energies nana im. 
In fact, HH and FY results for AV18+UIX model agree on cr(0) = 1.73 b total cross section value, which is consistent 
with experimental one <r(0) = 1.70 ± 0.03 b of Ref. j^. Three nucleon forces, and UIX in particular, have however 
very small effect in the resonance region as shown in Fig [2 it has even the tendency to diminish the cross section due 
to triton rescaling |3l l30l l31|. The very weak sensitivity of total cross sections to "standard" 3NF, may indicate that 
this disagreement could have its origin in the NN forces themselves and their P-waves in particular [lij . 

Let us remember that in N — d scattering there exist other large discrepancies between theory and experiment, the 
most known example is the so called A y puzzle Erl l47l |4S| . To explain these discrepancies, speculations about either 
the deficiency of the NN potentials in 3 Pj waves 49], or the presence of exotic 3N force terms not contemplated so 
far [50L l5l|| . have been advanced. Whether such discrepancies in 3N and 4N systems originate from the same deficiency 
in the nuclear Hamiltonian is still under debate. 

Some doubts could be casted on the reliability of the experimental data in the resonance region, due to existence 
of a single accurate data set |4^. However, due to the near agreement with calculations at low energies, its hardly 
believable that experimental points could contain a normalization error near the resonance peak. A normalization 
error can also be excluded by analyzing differential cross sections 01 : they indicate that theoretical results agree 
with experimental ones at cm = 90 deg, where positive parity phaseshifts dominate, while diverge at forward and 
backward directions, where negative parity phaseshifts become important. Therefore the experimental data of ref. 0] 
seems to be accurate and coherent. In addition, at near-resonance energies, this data seems to be also in agreement 
with older and less accurate measurements [H^ . Nevertheless, in view of the very few experimental studies of n- 3 H it 
would be very interesting to have an independent confirmation of the data. 
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FIG. 2: Comparison between experimental and theoretical n- H total cross section calculated using different methods and with 
AV18 potential. Dotted curve presents AV18+UIX model results of Ref. Q. 
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A similar analysis has been performed with FY and AGS methods using Nijm-II model. The obtained results are 
summarized in Table HVI Note that a particular choice of the local NN-interaction model has no qualitative impact. 
Thus Nijm-II predictions are very close to the AV18 ones and lead to the same discrepancies. 

TABLE IV: Nijm-II n — 3 H S- and P-wave phase-shifts (degrees) as a function of the center of mass kinetic energy E cm (MeV). 
Numbers in brackets are the mixing parameters in degrees for corresponding scattering states. 





0+ 


It 


1+ 


0" 


ir i 2 " 


T 


a(b) 




0.40 


-27.81 


-25.03 


-0.0183 


2.28 


2.84 3.89 


3.47 


1.73 


AGS 




-28.66 


-25.03 


-0.0167 (0.196) 


1.751 


2.633 3.537 (-39.44) 


3.349 


1.74 


FY 


0.75 


-37.27 


-33.68 


-0.072 


5.67 


6.63 9.94 


8.92 


1.79 


AGS 




-38.57 


-33.73 


-0.696 (0.354) 


4.252 


6.040 8.763 (-41.03) 


8.505 


1.77 


FY 


1.50 


-50.68 


-46.11 


-0.271 


14.66 


14.98 25.83 


23.87 


2.21 


AGS 




-52.48 


-46.22 


-0.334 (0.661) 


10.74 


13.20 21.72 (-42.69) 


22.37 


2.04 


FY 


2.625 


-63.75 


-58.38 


-0.639 


27.77 


24.51 44.66 


43.69 


2.51 


AGS 




-66.49 


-58.51 


-0.813 (1.067) 


20.22 


20.44 36.96 (-43.31) 


40.82 


2.23 


FY 


3.0 


-67.18 


-61.81 


-0.790 


32.00 


26.77 48.13 


48.04 


2.47 


AGS 




-70.04 


-62.67 


-1.10 (1.195) 


23.09 


22.13 40.35 (-43.29) 


45.25 


2.20 


FY 



IV. CONCLUSION 

In this work we have presented the low energy n — 3 H elastic cross sections calculated by solving the 4-nucleon 
problem with realistic NN interactions and by using three different approaches. 

The results of Hyperspherical Harmonics and Faddeev-Yakubovsky methods are converged to 1% accuracy and are 
in close agreement with each other. AV18 and Nijm-II nucleon-nucleon potentials underestimate the experimental 
cross section by 10% near the resonance peak. This disagreement is not corrected by UIX three-nucleon forces, whose 
effect in this energy region is small and tends to further decrease the theoretical value upto 15% |30l l3ll l44|. 

The results of AGS equations are closer to the experimental data. However they have been obtained using a 
rank-one approximation to the two-body t-matrix. This approximation, which was found to provide excellent results 
for nuclear systems dominated by S-waves (3N problem or positive parity 4N states), seems to provide an excess of 
attraction in the negative parity states of n — 3 H system, where P-waves turn out to be important. We believe that 
accurate treatment of these resonant states requires the inclusion of higher rank terms. 

The results presented here concern local NN potentials. However similar conclusions were found in previous 
works |3| using non local ones (Bonn |5^| and Doleschall [54|). This indicates a serious difficulty of the existing 
NN force models in describing the simplest nuclear resonance, i.e. the n — 3 H system. This difficulty can hardly be 
solved by the inclusion of a standard type of 3NF, used to reproduce the few-nucleon binding energies. Their origin 
could rather lie either in the NN forces themselves, or in the presence of 3NF of unknown type. 
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